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Heuristic Solving Linear Programming Problems 
Heurystyczne rozwiązywanie problemów programowania liniowego 


Эвристическое решение проблем линейного программирования 


The exact solutions of linear programming problems are obtained 
by means of simplex method and the electronic calculating machines. 
The paper is aimed at showing that finding the approximate solutions 
of certain important practically linear programming problems by means 
of heuristic approach requires using neither the sophisticated algorithm 
nor the electronic computer. Let us notice that the advisability of such 
an approach to solving linear programming problems as well as the 
practical usefulness of the approximate solutions had been first empha- 
sized by L. V. Kantorowicz (1). 

The paper is divided into two parts. The general exposition of the 
heuristic approach is followed by several examples of its application. 


1. Let us consider the following linear programming problem: 
maximize the objective function 


$ = XTR (1) 
subject to 
АХ-Ү--Р, (2) 
X20 (3) 
Y>0 (4) 
where 


X — the column vector (nX 1) of the decision variables, 
XT — the transposition of the vector X, 
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Y — the column vector (їп Х 1) of the slack variables, 

В — the column vector (nX1) of the objective function parameters, 
A — the matrix (п n) of the coefficients, 

Po — the column vector (mX 1) of the constants. 


Assuming that some coefficients in the matrix A are equal to zero! 
the approximate solution of the problem (1)—(4) can be found by means 
of the following heuristic procedure: 

a) Finding the first solution: 


— choosing the first sequence of the decision variables: 


` Kal KE ceert ын 
оп the basis of 
a ees > г} 
_ where 
]ЕМ 
N — the set of the decision variables subscripts, 


rj — the j-th element of the vector В. 


— finding the maximum value, of the variable x} : 


1.7 Н | b, b; $ b, | 
max, x; = тіп. |--, — ....... 
aij azj ал, 
\уһеге 
aij — the element of the matrix A standing in the i-th row and the 


j-th column, 
b; — the i-th element of the vector Ро. 


- finding the vector Р,, of the constants: 
Ры =P,—x} Pj 


where 
Р; — the j-th column vector of the matrix A corresponding with the 
variable xj | 


1 [et us notice that it das not pay to use the heuristic approach when there 
are all coefficients bigger than zero in the matrix A. - 


/ 
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- finding the maximum value of the variable x?: 
b; b; bs 
max. x? = ЕЗ "E m | 
aiy аз) Dms 
where 
b; — the i-th element of the vector Р, 
— finding the vector Pi, of the constants: 
Po, = Pox; Р? 
where 
P? — the j-th column vector of the matrix A corresponding to the 
variable x? 
— finding the maximum value of the variable x ? 
- b: b 
тах. x? = mtn (> 225 $E ыры m ) 
АХЛТ аз; amj 
where 
b;"" — the i-th element of the vector Р; “~. 
— finding the vector Pj" of the constants: 
р Р, өө PO vee —x? Р? 
where 
P? — the j-th column vector of the matrix A corresponding to the 
variable x? 


— finding the values of the slack variables: 
Y = р 


b) Finding the second solution: 


— choosing the second sequence of the decision variables on the ba- 


sis of the first sequence and the results of the first solution ? 


an example below. 


2 The way of performing that operation will be further explained when solving 
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— performance of the operations shown above when looking for the 
first solution. 


(s) Finding the s-th solution ? in the way shown above. 
(s+1) Choosing the approximate solution of the problem (1)—(4): 
— finding the approximate maximum value of the objective function 
on the basis of 
appr. max. S = max. (!XTR, ?XTR, ..., SXTR) 
where 
PXT (р = 1, 2,...,8) — the vector (1хп) of the values assigned 
to the decision variables in the p-th solution 
— defining the approximate solution: 


assuming that the appr.max.S has been achieved in the p-th solution 
then the approximate solution is given by the vectors: PXT and 
У = Po,” where Po,~ is the vector of the constants obtained in the 
p-th solution 3. 


An example. 


Maximize the objective function 


3 ч 
ј=1 
subject to 

811X1 + аззхз+ ajX4 T yi = by 
a22X2 + у» = be 
AgeXe + азах4 "E ya = b; 
841X4 + Уг = ba 
+ аѕзхз + 854Х4 + ys = bs 


xji20 (j=1..4) 
Yj+k 20 (k — Ts 9) 


3 As our experience shows, at least two solutions are necessary for choosing an 
approximate solution satisfactory practically. 
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(a) Finding the first solution: 


— assuming that 274 >> 23 > 21 Zə, then the first sequence of the 
decision variables is as follows: 


X4 Хз Ху Xe 
— assuming also that 


: | b, b, b; | 1 
max. x, = тіп. |——, —, — |= — 
ais аз 854 814 
70 
b, 814 ba b 
| be b, 0 bs—ag, — 
= then Ро, = Р|,-ХаР, == b, —— Aga = 214 
b, ам 10 b, 
b; 854 bi 
b,—ag4 — 
- aw — 
b, 
7 0 b,—a4, —— 
— max. хз = min. | — , ам = 0 
813 eo 
853 
mo =) 
— max. X, = min. | —, — |= 0 
811 841 
b b 
bs b,—a 4 БЕА ` b -а ZA 
Е : 3 з — аза 
— assuming that max. x; = тіп. | — , 814 | = 814 
22 m n 
0 
b, b, 0 
ы-ы | slm 
— then Pa = Ро —x;P, = "aa |- м [аз | = 
5, азг 0 
р 0 
bs—a54 — 
314 _ 
0 MI 
ba —às2X; Ya 
-10 — | Уз 
b, Ya 
bs —àg4X4 Ys 
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(b) Finding the second solution: 


— as the second sequence of the decision variables we assume: 
X3 X4 Xj Xo 


because in the first solution хә 2 0 and x= х; = 0 in spite of that 
73 С> 24 2 Zg. | 


— assuming that 


: | b, 5 ) b 
max. Xa = min.|—, — | = — 
813 853 453 
b, 813 b, —213X3 
b, b 0 b, 
— then Po, —PxP,—]|bs|— 5 0 = | b; 
b, | 2910 b, 
b; а5з 0 
: | bı —a13X3 bs 0 | 
— max. ха = min. |, ==, —|=0 
а14 334 a54 
— assuming next that 
2 : | b, —a45X3 b, | 4 
max. x, = min. {--—— --|--- 
814 ад ад 
6,—аџз X3 811 bı—(a13X3 +a11x1) 
, "n b, b, |? b, 
— thenP,, = P,,—x,P, | bs —— 0 = | bə 
b, ап ац 0 
0 0 0 


— and assuming lastly that 


| а be 
тах. x3 = min. |—— —] = 


Ane 333 аз 
b, — (а1зхэ-Назах,) 0 
5i 3 1 Ь, даа 
— then Poo = Ро, —х,Р, = b, = or Ago = 
0 2210 
0 
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b; — (а1зхз-Раха) y 
0 Ye 
= | Бз—азох» —] Уз 
0 | Ул 
0 ys 


(c) Choosing the approximate solution: 


— assuming that 
appr. max. S = max. (IXTZ, 2XTZ) = ?XTZ 
— then the approximate solution is as follows: 


| b, b; . bs 
Хрэ--  Хос--- OH = — yy = В: — (а1зхз-Рал 1х1) Ya = Бэ--азаХа 
ад 853 85g 


Let us notice that there are exactly five values larger than zero in 
the approximate solution which signify that the solution is basic. 


2. We shall now turn to the applications of the presented heuristic 
approach. The nature of the first problem we are going to deal with is 
as follows: 


to find such production program XT = [X;, хо, ..., Xy] which maximizes 
approximately the objective function 


E =-= 
X ny, 

i =l 

and fulfills the conditions 
AX «Р, 
X20 
where 
му, пу — the unit output, input, respectively 


ү 

J : А 
— — the unit effectiveness 
n 

J 


19* 
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w 
= — — the average effectiveness 
n 


where 


w, n — the average output, input, respectively 


Thus 
wY К w,Y 
—|zE*'zl—- 
n, n, 
where Ue m c 


Еа — the maximum value of the objective function, 
w,\ /w,\” | 
ксы Wee ae the smallest, largest unit effectiveness, respectively 
1 1 
Let us now consider the following numerical example: . 
maximize the objective function 
37x, 4-28x4 4- 30x; 4-25x, 4- 17x; 4- 18x, 
77 148x,+122x9+140x 3+ 136x4+ 100x; +120x, 


- - 


subject to 
Xi у: = 8 000 
хә +ye = 25 000 
90x, + 80x2 + 75x3+ 50x4+ 52x5+ 60x6 +уз = 2 870 000 
30x1 + 20х + 20ж- 18x44- 46х5+ 34xq +у, = 775000 


xj (j= 1,2...6) 20 
Ук(К--1,2..4) 20 


(a) Finding the first solution: 
— the first sequence of the decision variables is 


Ху Xo Хз X4 X5 Xe 
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because the unit effectiveness indicators г; are 
j 12 3 4 5 6 
г) 25 23 21 18 17 15 
where 


гу = [wj/nj] 100 


RA | (296 28700008) 1775000 
— the . X, = min. > ~o ~ ; Е 
п тах. х; m 1 90 30 
— 25 833) — 8000 
8 000 | и | 0 
‚ 25 000 0 25000 
— then Py = Р.м, = | 5 870000 | —8000 190175 | 2 150000 
775 000 30 535000 
25000 2150000 535000 
ыг A = in. , —.— = 26875 ? 70 — 
then max. xs — min | 1 80 20 
= 26750) = 25000 
: х | j 
р 250 00 1 0 
— then B, = Po,—x,P, = 2 150 000 —25000 80| — 1150000 
535 000 20 ‚ 35000 
150000 35000 
— then max. хз — min. | == 20000, 202 1750) = 1750 
0 0 0 
MC 0 0 0 
— th = Py, —%;P; = = = 
en Pa = Р„—х,Р, 150000 1750 75 18 750 
_ 35000 20 0 
y 0 
= ст Уг. 0 
then x, = № 67-0 and | — [ 615 
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As can be inferred from the foregoing discussion the solution. just 
obtained can be already considered as the one looked for. The corres- 
ponding value of the -objective function is the maximum value Ea, 
assuming that the values of the slack variables are minimized. 

The next problem to be considered here is given in Table 1. When 
solving it, let us first notice and take advantage of its specific structure. 
Namely, the set of 16 decision variables can be divided into two parts. 
To the first part belong those variables values of which may be deter- 
mined independently on the basis of only one appropriate constraint: 


х; . (2) 
Xe, X7, Xs, Xio, X11, X16 (3) 9 
X12 (5) 


and to the second part those remaining, values of which cannot be found 
in such a way: 


(a) Finding the values of the decision variables belonging to the first 
category: 


120 ; 
— тах. ху = ==2 and уз=0 


322 


— max. Xi = 322 = 


10 and y;=0 


— the sequence of the variables xe, хт, Xs, X10, Хи, Х is as follows: 
§ 16 6 7 8 11 10 
because the standardized parameters of the objective function are 


j 6 7 8 10 11 16 
ey 0,028 0,025. 0,0224 0,016 0,0215 0,085 


where e, = гра»: 


where г; — the j-th parameter in the objective function, 
азу — the j-th coefficient in the second constraint. 


2430 
— hence max. Xie = ET = 


30 


and xg = x, = xg = X49 = хи = 0 


and у; = 0 


4 The reference number of the constraint, 
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(b) Determining the values of the variables belonging to the second 
category: 


— to maximize 
Е’ == 0,88х, + 1,94х,+ 3,8 хз + 5,44х,-+ 6,86 хо + 1,94х;3+ 10,5хџ + 10,93х;;5 
subject to 


(1) 89,3x,+ 164х,+ 132,2x4 + 199,7 x4-- 501х»-+ 67,6х1з-+ 603,6x14 + 469,3x15 + 


(4) Б 49, 93x, + 26,04х, + 26,04х;5+ 
(6) 16,11x + 33,83x, 
(7) 16x,4- 

+ Ут == 20 


— the sequence of the decision variables is 
Хз X13 X4 Хуб Xia Хэ Хә Xi 
because the standardized parameters of the objective function аге 


j 1 2 3 4 9 13 14 15 
f, 0,010 0,012 0,03 0,027 0,014 0,028 0,017 0,023 


where 1, = гу/а; 


where ау; — the j-th coefficient in the first constraint 


— mi | 221120 =217) = 217 
атта 
3 424 132,2 3 137,13. 
46 0 46 
— = -8=| МТ ent | о 


20 0 | 20 
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(0772 1= ҹу) оёх 


Í (91: zt=9 оз 
02 = f+ 41х91 
ce + “a *xegtes+ "xr ror 
we = А+ i nyg'ze 
gp = S xy0'9z -- "xv0'97 Ч +-°xe6'6r 
021 = + +%х09 
ogre = FA + Өх184-"х9014-55Х8184- +PXT'ZE+iX6 ZZ++’ ЕТ 
ҮРЕ =M E? x£'695 4- x9'c09 4-*1x9'19 1- J-*x toc -'X1'66T--"Xz'c£T--*xv9T--*xe'68 


03 yoefqns 
*Tx16'9 4-3 xe6'01 - xc'01 + xpe'T 4-5 х8с--7х8024-01х6"6--0х98 9--5Х21/0--2Х86/0 --*x8y'0 --?xcc'g 4-"xpy'c ХВЕ --*xv6' T--"xgg'o = Я 


uonounj eAnjefqo ay} eziurixeul OL 


“савд, 


(9) 
(8) 
(9 
(9) 
(9) 
(p) 
(=) 
(2) 
(D 
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2 3137,13 odi 
== тах. X13 = 67,6 == , 
3137,13 675 0 
: ; 46 0 46 
- р = Р-Р | y |-404161-10 
20 o | [2 


— hence x, = хо = X4 = Ху = xq, = X15 = 0 


у 1 0 

Ya 46 
and = 

Ye 0 

Уг 20 


Taking into account the shown specific structure of the problem апа 
the performed standardization of the objective function parameters we 
shall assume as the approximate solution the one just obtained: 


хз = 2,17 X; = 2 X12 = 10 X43 — 46,41 Xig = 30 


у: = 46 у: = 20 


Let us notice moreover that the outlined here general idea of the 
heuristic solving linear programming problems has — as our experience 
proves — many applications of practical importance, e. g. in finding the 
production assortment at the factory level. 
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STRESZCZENIE 


_Celem tego artykulu jest ukazanie, ze aby znalezé przyblizone rozwiazania pew- 
nych praktycznie istotnych problemów programowania liniowego та pomoca po- 
dejscia heurystycznego, nie jest konieczne stosowanie ani skomplikowanego algo- 
rytmu, ani komputera elektronicznego. 

Artykuł dzieli sie na dwie części. Po ogólnym przedstawieniu proponowanej 
procedury heurystycznej nastepuje kilka przykładów jej zastosowania, 


РЕЗЮМЕ 


Цель статьи — показать, что для того, чтобы найти приближенпое решение 
некоторых существенных с практической точки зрения преблем минейного про- 
граммирования с помощью эвристического подхода, необязательно применять ни 
сложный логарифм, ни электронный компьютер. 

Статья состоит из двух частей. После общего представления предлагаемой 
эвристической процедуры приводятся примеры ее применения. 


